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We analyse a new mechanism for the cosmological QCD first-order phase transition: 
inhomogeneous nucleation. The primordial temperature fluctuations are larger 
than the tiny temperature interval, in which bubbles would form in the standard 
picture of homogeneous nucleation. Thus the bubbles nucleate at cold spots. We 
find the typical distance between bubble centers to be a few meters. This exceeds 
the estimates from homogeneous nucleation by two orders of magnitude. The 
resulting baryon inhomogeneities may affect primordial nucleosynthesis. 



a To appear in the Proceedings of Strong and Electroweak Matter 2000 (SEWM2000), 
Marseilles, France, 14-17 June 2000. 



First-order phase transitions normally proceed via nucleation of bubbles of 
the new phase. When the temperature is spatially uniform and no impurities 
are present, the mechanism is homogeneous nucleation. We denote by Wheat the 
effective speed by which released latent heat propagates in sufficient amounts 
to shut down nucleations. The mean distance between nucleation centers, 
measured immediately after the transition completed, is 

^nuc,hom — 2l>h ca t At nuc , (f) 

where Ai nuc is the duration of the nucleation periodEH. 

Surface tension and latent heat are provided by lattice simulations with 
quenched QCD only, giving the values a = 0.0f5T c 3 , I = 1.4T*H. Scaling the 
latent heat for the physical QCD leads us to take I = 3T 4 . Using these values 
the amount of supercooling is A sc = 2.3 x 10 -4 , with large error bars. Due 
to rapid change of energy density near T c , the microscopic sound speed in the 
quark phase has a small value, 3c 2 (Tc) < 1 (for references seelil). 

In the real Universe the local temperature of the radiation fluid fluctuates. 
The temperature contrast is denoted by A = 5T/T. On subhorizon scales in 
the radiation dominated epoch, each Fourier coefficient A(f, k) oscillates with 
constant amplitude Ay(fc). Inflation predicts a Gaussian distribution, 

1 / 1 A 2 \ 
p(A)dA = -= exp — — ^ dA . (2) 

{ ' V2^A™ S V 2(A™ S )V 1 ; 

We find! for the COBE normalized! rms temperature fluctuation of the radia- 
tion fluid (not of cold dark matter) A™ s = 1.0 x 10~ 4 for a primordial Harrison- 
Zel'dovich spectrum. The change of the equation of state prior to the QCD 
transition modifies the temperature-energy density relation, A = c 2 fe/(e + p). 
We may neglect the pressure p near the critical temperature since p <C £ q at 
T c . On the other hand the drop of the sound speed enhances the amplitude 
of the density fluctuations proportional to c s i. Putting all those effects 
together and allowing for a tilt in the power spectrum, the COBE normalized 
rms temperature fluctuation reads 

/ L. \ (»-l)/2 

A™ s » 10" 4 (3c 2 ) 3 / 4 (jL j . (3) 

For a Harrison-Zel'dovich spectrum and 3c 2 = 0.1, we find A™ s w 2 x 10~ 5 . 

A small scale cut-off in the spectrum of primordial temperature fluctua- 
tions comes from collisional damping by neutrinosQ. For Z SBl ™>th) the length 
over which temperature is constant, we take the value 10~ 4 g?h^- The temper- 
ature fluctuations are frozen with respect to the time scale of nucleations. 
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Figure 1: Sketch of a first-order QCD transition in the inhomogeneous Universe. At t± 
the first hadronic bubbles (H) nucleate at the coldest spots (light gray), while most of the 
Universe remains in the quark phase (Q). At the bubbles inside the cold spots have merged 
and have grown to bubbles as large as the temperature fluctuation scale. At t% the transition 
is almost finished. The last quark droplets are found in the hottest spots (dark gray). 

Let us now investigate bubble nucleation in a Universe with spatially inho- 
mogeneous temperature distribution. Nucleation effectively takes place while 
the temperature drops by the tiny amount A nuc . To determine the mechanism 
of nucleation, we compare A lmc with the rms temperature fluctuation A™ ls : 

f . If A™ ls < A nuc , the probability to nucleate a bubble at a given time is 
homogeneous in space. This is the case of homogeneous nucleation. 

2. If A™ ls > A nuc , the probability to nucleate a bubble at a given time is 
inhomogeneous in space. We call this inhomogeneous nucleation. 

The quenched lattice QCD data and a COBE normalized flat spectrum 
lead to the values A nuc ~ 10~ 6 and A™ — 10~ 5 . We conclude that the 
cosmological QCD transition may proceed via inhomogeneous nucleation. A 
sketch of inhomogeneous nucleation is shown in Fig. [[} The basic idea is 
that temperature inhomogeneities determine the location of bubble nucleation. 
Bubbles nucleate first in the cold regions. 

For the fastest fluctuations, with angular frequency c s /Z smoo th, we find 



The Hubble expansion is the dominant contribution, as typical values are 3c^ = 
0.1 from quenched lattice QCD and A™ s £h/<^ ~ 0.01. This means that the 
local temperature does never increase, except by the released latent heat during 
bubble growth. 

To gain some insight in the physics of inhomogeneous nucleation, let us first 
inspect a simplified case. We have some randomly distributed cold spheres of 




diameter Z sm0 oth with equal and uniform temperature, which is by the amount 
A™ S T C smaller than the again uniform temperature in the rest of the Universe. 
When the temperature in the cold spots has dropped to Tf, homogeneous 
nucleation takes place in them. Due to the Hubble expansion the rest of the 
Universe would need the time Ai coo i = i#A™ s /3Cg to cool down to Tf. The 
cold spots have fully been transformed into the hadron phase while the rest of 
the Universe still is in the quark phase. The latent heat released in a cold spot 
propagates in all directions, which provides the length scale 

^heat = 2fhcat At coo i. (5) 

If the typical distance from the boundary of a cold spot to the boundary of a 
neighboring cold spot is less than Zheatj then no hadronic bubbles can nucleate 
in the intervening space. 

The real Universe consists of smooth patches of typical linear size Smooth, 
their temperatures given by the distribution (j^). At time t heat, coming from 
a cold spot which was transformed into hadron phase at time t 1 , occupies the 
volume V(t, t 1 ) = (47r/3)[/ smoo th/2 + fhcat(i — t')] 3 - The fraction of space that 
is not reheated by the released latent heat (and not transformed to hadron 
phase), is given at time t by 

/(<)«! - [ T ihn (t')V(t,t')dt', (6) 
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where we neglect overlap and merging of heat fronts. T^n is the volume fraction 
converted into the new phase, per physical time and volume as a function of 
the mean temperature T = T(t). We find 

r ih „ = 3c^— - p(A = ^ - 1), (7) 

J EH ^smooth J 

where the relevant physical volume is V smoo th = (47r/3)(7 smooth /2) 3 . 

The end of the nucleation period, i;hm is defined through the condition 
f(t ihn ) = 0. We introduce the variables N = (1 - T f /T)/A™ s and Af = 
N(tih n ). Putting everything together we determine Af from 

JL^ r dN ?^L {b2Z2* +N _jj\ =1 (8) 
'smooth J N V27T V 'heat / 

For Zhcat/'smooth = 1, 2, 5, 10 we find Af « 0.8, 1.4, 2.1, 2.6, respectively. 

The effective nucleation distance in inhomogeneous nucleation is defined 
from the number density of those cold spots that acted as nucleation centers, 
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dtmcihi, = n 1/3 . We find 



^nuc.ihn 



T ihn (t)dt 



i) 



-1/3 



(9) 



= [-(1 - erf(AA/\/2)]- 1/3 'smooth- (10) 

With the above values Zhcat/'smooth = 1,2,5, 10 we get d nuC!ihn = 1.4, 1.8,3.0, 
4.8 x Smooth, where Smooth ~ 1 m. 

For a COBE normalized spectrum without any tilt and with a tilt of 
n — 1 = 0.2, together with 3c^ = 0.1 and Wheat = 0.1, we find the estimate 
'heat /Smooth ~ 0.4 and 9, correspondingly. Notice that the values of Wheat 
and ic 2 s are in principle unknown. Anyway, we can conclude that the case 
'heat > 'smooth is a realistic possibility. 

With 2wheat(3Cg) -1 / 4 (10 -4 (iH/'smooth) < 1 and without positive tilt we are 
in the region /heat < 'smooth, where the geometry is more complicated and the 
above quantitative analysis does not apply. In this situation nucleations take 
place in the most common cold spots (J\f ~ 1), which are very close to each 
other. We expect a structure of interconnected baryon-deplcted and baryon- 
enriched layers with typical surface 's moot h an d thickness Zdof = Wdcf A^cooi. 

We emphasize that inhomogeneous and heterogeneous nucleationQ are gen- 
uinely different mechanisms, although they give the same typical scale of a few 
meters by chance. If latent heat and surface tension of QCD would turn out 
to reduce A sc to, e.g., 10 -6 , instead of 10 -4 , the maximal heterogeneous nu- 
cleation distance would fall to the centimeter scale, whereas on the distance in 
inhomogeneous nucleation this would have no effect. 

We have shown that inhomogeneous nucleation during the QCD transi- 
tion can give rise to an inhomogeneity scale exceeding the proton diffusion 
scale (2 m at 150 MeV). The resulting baryon inhomogeneities could provide 
inhomogeneous initial conditions for nucleosynthesis. Observable deviations 
from the element abundances predicted by homogeneous nucleosynthesis seem 
to be possible in that case a. 

In conclusion, we found that inhomogeneous nucleation leads to nucleation 
distances that exceed by two orders of magnitude estimates based on homo- 
geneous nucleation. We point out that this new effect appears for the (today) 
most probable range of cosmological and QCD parameters. 
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